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Abstract - We study the superfluidity of a spin-orbit coupled Bose-Einstein condensate (BEC) by 
computing its Bogoliubov excitations, which are found to consist of two branches: one is gapless 
and phonon-like at long wavelength; the other is typically gapped. These excitations imply a 
superfluidity that has two new features: (i) due to the absence of the Galilean invariance, one can 
no longer define the critical velocity of superfluidity independent of the reference frame; (ii) the 
superfluidity depends not only on whether the speed of the BEC exceeds a critical value, but also 
on cross helicity that is defined as the direction of the cross product of the spin and the kinetic 
momentum of the BEC. 



Superfluidity was first discovered in the 1930s, and has 
fascinated physicists ever since. This interesting phe- 
nomenon was explained by Landau [1] , whose theory has 
been very successful in explaining many important prop- 
erties of superfluids. However, Landau's theory of super- 
fluidity may be facing challenges brought by the recent 
experimental realization of artificial gauge fields for ultra- 
cold bosonic atoms |^-[5]. When the artificial gauge field 
is non-Abelian [BHH], it is effectively spin-orbit coupling 
(SOC). SOC has played a crucial role in many exotic phe- 
nomena such as topological insulators [9]. However, in su- 
perfluids, the SOC is generally absent and its effects have 
remained largely unexplored. Note that this issue is not 
limited to ultracold atoms since the Bose-Einstein conden- 
sation with SOC also exists for excitons in semiconductors 

mm- 

There have been some theoretical works, where many 
interesting properties of spin-orbit coupled BECs are ex- 
plored For instance, it was shown in Ref. [TT] that 
SOC can lead to unconventional Bose-Einstein condensa- 
tion with the breaking of time-reversal symmetry. Later 
a stripe phase that breaks rotational symmetry was found 
P^[T^ . In this Letter we concentrate on the superfluidity 
of the spin-orbit coupled BEC. 

To put our study into perspective, we briefly review 
Landau's theory for a superfluid without SOC. Consider 



such a supcrfiuid fiowing in a tube. With the Galilean 
transformation. Landau found that the excitation of this 
flowing superfluid is related to the excitation of a motion- 
less superfiuid as [T| 



e„(p) =eo(p) + p-v, 



(1) 



where p is the momentum of the excitation. For phonon 
excitation eo(p) = c|p|, the excitation e„(p) can be nega- 
tive only when |v| > c. Therefore, the speed of sound c is 
the critical velocity beyond which the flowing supcrfiuid 
loses its superfluidity and suffers viscosity. We switch to 
a different reference frame, where the superfluid is at rest 
while the tube is moving. It is apparent to many that these 
two reference frames are equivalent so that the superfluid 
will be dragged along only when the tube speed exceeds 
the speed of sound c. However, this equivalence is based on 
that the superfluid is invariant under the Galilean trans- 
formation. As SOC breaks the Galilean invariance of the 
system [24j , we find that these two reference frames are no 
longer equivalent as shown in Fig. [T] the critical speed for 
(a) is different from the one for (b). For easy reference, the 
critical speed for (a) is hereafter called the critical flowing 
speed and the one for (b) the critical dragging speed. 

Our study of the superfluidity for a BEC with SOC is 
based on the computation of the elementary excitations 
using the Bogoliubov equation. We calculate how the ele- 
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(b) 



SOC superfluid 




(c) 



Fig. 1: (a) A superfluid with spin orbit-coupling moves while 
the tube is at rest, (b) The superfluid is dragged by a tube 
moving at the speed of v. (c) An impurity moves at v in the 
SOC superfluid. The reference frame is the lab. The two-way 
arrow indicates the equivalence between different scenarios and 
the arrow with a bar indicates the non-equivalence. 



mentary excitations change with the flow speed and man- 
age to derive from these excitations the critical speeds for 
the two different scenarios shown in Fig. [TJa,b). We find 
that there are two branches of elementary excitations for 
a BEC with SOC: the lower branch is phonon-like at long 
wavelengths and the upper branch is generally gapped. 
Careful analysis of these excitations indicates that a mov- 
ing BEC with the SOC (Fig. [TJa)) has non-zero critical 
velocities for superfluidity while the critical speed is zero 
for a BEC being dragged by a moving tube(Fig. [Ijb)). 
This shows that the definition of critical velocity depends 
on the reference frame for a superfluid with SOC and, 
probably, for any superfluid that has no Galilean invari- 
ance. 

In addition, we find that the properties of a flow of 
BEC with SOC are also related to its spin direction. We 
characterize this spin direction with cross helicity, which 
is defined as the cross product of the spin and the kinetic 
momentum of the flow. A BEC flow with Rashba SOC is 
always unstable if its cross helicity is negative. 

Model: We consider a BEC with pseudospin 1/2 and 
the Rashba SOC. The dynamics of the system can be de- 
scribed by the Hamiltonian P^[T51[^[^ 
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^(|*l|4 + |vI/2|4)+C2|*invI/2p^ (2) 



Ci 
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where 7 is the SOC constant, Ci and C2 are interac- 
tion strengths between the same and different pseudospin 
states, respectively. 

We focus on the homogeneous case V{r) = despite 
that the BEC usually resides in a harmonic trap in exper- 



iments. The primary reason is that the superfluidity can 
be discussed more clearly in the homogeneous case, and 
be compared directly with the conventional superfluidity 
of a spinless bosonic system. In addition, the results in 
the homogeneous ease can be adopted to understand the 
superfluidity in more complicated situations with the lo- 
cal density approximation. We also limit ourselves to the 
case Ci > C2, where the system is stable against phase 
separation [T^[TS]. In the following discussion, for sim- 
plicity, we set h = m = 1 and ignore the non-essential z 
direction, treating the system as two-dimensional. This 
does not impair the validity of our model, since in the fol- 
lowing we assume the BEC moves in the y direction, and 
the critical velocity is found to be not influenced by the 
excitation in the z direction. 

The Gross-Pitaevskii equation obtained from the Hamil- 
tonian ([2]) has plane wave solutions 
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where tan6ik = k^/ky, ^(k) |k| 2/2 - 7|k| + (Ci + C2)/2. 
The solution 0k is the ground state of the system when 
|k| — 7. There is another set of plane wave solutions, 
which have higher energies and are not relevant to our 
discussion. 

The plane wave solution 0k represents a BEC flow with 
the velocity v = k — 7k(k = k/|k|). This velocity is essen- 
tially the kinetic momentum of the BEC, which is different 
from the conjugate momentum k due to the presence of 
SOC. For example, the ground state has a non-zero conju- 
gate momentum k. but its kinetic momentum (or velocity) 
is zero. Besides its velocity, the BEC flow described by (/)k 
has another feature, the direction of the spin, which has to 
be included for a complete description of the flow. For ex- 
ample, the flow at k = 87/2^ has the same velocity 7/2^ as 
the flow at k = —"f/2y. However, they have different spin 
directions. With the Rashba type SOC, the spin direction 
of the eigenstate is always perpendicular to the velocity 
and has only two choices. As a result, it is sufficient and 
also very convenient to use w = sign[(v x a) ■ z] to de- 
note the spin direction. This variable w, which is either 1 
or —1, is the cross helicity mentioned in the introduction. 
Note that the usual helicity in literature is defined as the 
inner product of the spin and the momentum. 

Critical velocities: We study first the scenario depicted 
in Fig. [TJa), where the BEC fiows with a given velocity. 
Since the system is not invariant under the Galilean trans- 
formation, we cannot use Eq. ^ to find the excitations 
for the flowing BEC from the excitation of a stationary 
BEC. We have to compute the excitations directly. This 
can be done by computing the elementary excitations of 
the state 0k with the Bogoliubov equation for different 
values of k. 

Without loss of generality, we choose k = /cy with k > 
0. Following the standard procedure of linearizing the 
Gross-Pitaevskii equation |27 1 f28]. we have the following 
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Fig. 2: (color online) Elementary excitations of a BEC flow with the SOC. (a) k = 1; (h) k — 2.5; (c) k — A. Ci = 10, C2 — 4, 
7 = 1. 
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These results show that the system at the ground state 
(fc = 0) has two different speeds of sound, ^ (Ci + C2)/2 
and -^Z (Ci — C2)/2. Since the excitation £_ becomes neg- 
ative only when k > \J [C\ — 02)/'^, the critical velocity 
in this case is ^ (Ci — C2)/2. When C2 = 0, these two 
branches of excitations merge into one and the critical ve- 
locity is yX?i/2, which is well known and was confirmed 
in a BEC experiment [30] . 

In general, there are no simple analytical results. We 
have numerically diagonalized C to obtain the elementary 
excitations. We find that part of the excitations are imag- 
inary for BEC flows with |k| < 7. This means that all the 
flows with |k| < 7 are dynamically unstable and therefore 



do not have superfluidity. For other flows with fc > 7, 
the excitations are always real and they are plotted in 
Fig. [21 One immediately notices that the excitations have 
two branches, which contact each other at a single point. 
More closer examination shows that the upper branch is 
gapped in most of the cases while the lower branch has 
phonon-like spectrum at large wavelength. These features 
are more apparent in Fig. [31 where only the excitations 
along the x direction and y direction are plotted. 



liiqx - qy + k) ~ and 643 = -l{iqx + qy - k) ~ H^. 
As usual, there are two groups of eigenvalues and only ^ 
the ones whose corresponding eigenvectors satisfy — 
\vi\^ = 1 (i = 1, 2) are physical. 

For comparison, we consider the case without SOC. This 
is to put 7 = and reduce our system to a two-component 
BEC system, which is well studied in literature [ISj- For 
this case, C can be diagonalized analytically and there are 
two branches of excitations 




Fig. 3: Excitations along the x direction (the first row) and 
y direction (the second row) at different values of k. (al, a2) 
= 1; (bl, b2) = 3; (cl, c2) fc = 4. Ci = 10, C2 = 4, 7 = 1. 

In Fig. [^Jc) and Fig. [31^c2), we notice that part of 
the excitations in the upper branch are negative, indicat- 
ing that the underlying BEC flow is thermodynamically 
unstable and has no superfluidity. In fact, our numeri- 
cal computation shows that there exists a critical value 
kc- when k > kc either part of the upper branch of excita- 
tions or part of the lower branch or both become negative. 
This means that the flows described by the plane wave so- 
lution (/)k,- with |k| > kc suffer the Landau instability 
and have no superfluidity. Combined with the fact that 
the flows with |k| < 7 are dynamically unstable, we can 
conclude that only the flows with 7 < |k| < kc have super- 
fluidity. Physically, these super-flows have speeds smaller 
than Vc = kc — "/ and cross helicity w = I. We have plot- 
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Fig. 4: Critical flowing velocity Vc (solid line) of a BEC as a 
function of the SOC parameter 7. (a) Ci > 3C2; (b) Ci < 3C2. 



ted how the critical flowing velocities vary with the SOC 
parameter 7 in Fig. 21 

We find two different asymptotic behaviors for the crit- 
ical velocity when 7 gets large. As shown in Fig. |4l 
the critical flowing velocity becomes a constant, Vc — 
\J (Ci + C2)/2, beyond a threshold value in the param- 
eter regime C\ > SC2 while it approaches an asymptotic 
value \^'TJi'^^C2 in the regime Ci < 36*2 . 

We turn to another reference frame illustrated in Fig. 
[2b), where the BEC can be viewed as being dragged by 
a moving tube. To simplify the discussion, we replace the 
moving tube with a macroscopic impurity moving inside 
the BEC as shown in Fig. [IJc) . Correspondingly, the ques- 
tion "whether the BEC will be dragged along by the mov- 
ing tube?" is replaced by an equivalent question "whether 
the impurity will experience any viscosity?" . Suppose the 
moving impurity generates an excitation in the BEC. Ac- 
cording to the conservations of both momentum and en- 
ergy, we should have 



moVi 



mov/ 



+ eo(q). 



(7) 
(8) 



where toq is the mass of the impurity, and Vf are the 
initial and final velocities of the impurity, respectively, and 
eo(q) is the excitation of the BEC at fc = 7. If the exci- 
tations were purely phonons, i.e., eo(q) = c|q|, these two 
conservations would not be satisfied simultaneously when 
u « |v.j| « |v/| < c. This means that the impurity could 
not generate phonons in the superfluid and would not ex- 
perience any viscosity when its speed were smaller than 
the sound speed. This is in fact nothing but the Cerenkov 
radiation where a charged particle radiates only 

when its speed exceeds the speed of light in the medium. 
For our BEC system, the elementary excitations eo(q) are 
not purely phonons. In this case, the critical dragging 
velocity derived from Eqs. (|7|8|) is given by 



eo(q) 



For the special case 7 = 0, we have 



ef(q) 
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(10) 



From Eq. (jH]), we obtain the critical dragging speed 
^^^2^^ , which is the same as the critical flowing speed. 
When 7 7^ 0, The excitations £o(q) also share two 
branches. Along the x direction, these two branches are 



(11) 

where si = + 7^ {Ci - C2), S2 = + \C^, h = s?, 

t2 = 2siS2, and ^3 = 2si + (7^ C2/2)^ Along the y 
direction, the excitations of the ground state are 



£o(%) = V ^ ^y + f' 



(12) 



4(g,) = 27g, + ^/2si+(.2-^)g2 + |.(13) 



When 7 > 0, the upper branch e^((ix) is always parabolic 
at small with a gap \J2s\. When expanded to the 
second order of qx-, the lower branch has the following 
form 



£0 (fe) 
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This shows that s^{(ix) is parabolic at long wavelengths 
instead of linear as usually expected for a boson system. 
This agrees with the results in Ref . [T7] . This parabolic ex- 
citation has a far-reaching consequence: according to Eq. 
(jni), the critical dragging velocity is zero, very different 
from the critical velocity for a BEC moving in a tube. This 
shows that the critical velocity for a BEC with SOC is not 
independent of the reference frame, in stark contrast with 
a homogeneous superfluid without SOC. This surprising 
finding of course has the root in the fact that the BEC de- 
scribed by the SOC Hamiltonian ([2]) is not invariant under 
the Galilean transformation [24] . 

Rashba and Dresselhaus SOC: We have also investigated 
the superfluidity with the general form of SOC, which is a 
mixture of Rashba and Dresselhaus coupling. Mathemat- 
ically, this SOC term has the form aaxPy — fSayPx- The 
essential physics is the same: the critical flowing speed is 
different from the critical dragging speed, and therefore 
the critical velocity depends on the choice of the reference 
frame. However, the details do differ when a ^ 13. First, 
the critical dragging speed is no longer zero. Secondly, 
the BEC with negative cross helicity can also be stable; as 
a result, there is a different critical speed for either cross 
helicity. For the Rashba type (a = /3), the critical speed 
for negative cross helicity is always zero. 

Experimental observation: Spin-orbit coupled BECs 
have been realized recently by three different groups [SHS] 
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through coupling ultra-cold ^^Rb atoms with laser fields. 
The strength of the SOC in the experiments can be tuned 
by changing the directions of the lasers [3H5] or through 
the fast modulation of the laser intensities The in- 

teraction between atoms can be adjusted by varying the 
confinement potential, the atom number or through the 
Fcshbach resonance For the scenario in Fig. [TJb), 

one can use a blue-detuned laser to mimic the impurity 
for the measurement of the critical dragging speed similar 
to the experiment reported in Ref. |30j . For the scenario in 
Fig. mja), there are two possible experimental setups for 
measuring the critical flowing speed. In the first one, one 
generates a dipole oscillation similar to the experiment in 
Ref. [5] but with a blue detuned laser in the middle of the 
trap. The second one is more complicated: At first, one 
can generate a moving BEC with a gravitomagnetic trap 
[55] . Afterward, one can use Brag spectroscopy [5^157] to 
measure the excitations of the moving BEC, from which 
the superfluidity can be inferred. For the typical atomic 
density of 10^* ~ 10^^ cm~^ achievable in current experi- 
ments [30], the critical flowing velocity is 0.2 ~ 0.6 mm/s, 
while the critical dragging velocity is still very small, about 
10~^ ~ 10~^ mm/s. To further magnify the difference be- 
tween the two critical velocities, one can use the Feshbach 
resonance to tune the s-wave scattering length. 

Perspective: With BEC, we now have many types of su- 
perfluids, which are different from the conventional spin- 
less and homogeneous supcrfluid typifled by liquid helium. 
These new superfiuids can be classified into two groups 
according to how much they differ from the conventional 
superfluid. We call the first group Landau superfiuid, for 
which Landau's argument (Eq. [1]) can be used to find 
its critical velocity; we call the second group non-Landau 
superfluid, for which Landau's argument is either insuffi- 
cient or invalid for finding its critical velocity. We offer 
two examples. 

The first is a two-component BEC without SOC. As 
discussed previously, although this superfluid has two 
branches of phonon excitations and appears very differ- 
ently from the conventional superfluid, its critical velocity 
can still be derived from the excitation at the ground state 
with Landau's argument. Therefore, this is a Landau su- 
perfluid. 

The second is a BEC in an optical lattice or supersolid 
helium. These periodic superfluids differ from the conven- 
tional superfluid by having a periodically modulated den- 
sity. We consider first the case when the BEC is locked 
with the lattice and they move together as shown in Fig. 
[5Ua,b,d). Since this system is invariant under the Galilean 
transformation, the two scenarios depicted in Fig. [S][a,b) 
are equivalent and Landau's argument is still applicable. 
The caveat is that the critical velocity is always zero no 
matter how the elementary excitation of the superfluid 
looks. The key reason is that the momentum p of the 
excitation is not well defined due to the presence of the 
lattice: two momenta differ by a reciprocal lattice vector 
are equivalent. This result has been verified in extensive 
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Fig. 5: In (a,b,d), tlie periodic superfluid and the lattice are 
locked in and move together relative to the obstacle. In (c), the 
lattice is gradually accelerated to the speed of v. The reference 
frame is the lab. The scenario in (c) is not equivalent to the 
scenarios in (a,b,d). 



computations for Cerenkov radiation in a periodic media 

The addition of a lattice does put a new twist into the 
system. One can gradually accelerate the lattice to a cer- 
tain velocity and see how the superfluid changes. There 
exists a critical velocity of the lattice beyond which the 
system loses its superfluidity. This critical velocity can 
no longer be determined by Landau's argument as having 
been discussed in Ref. [35]. Therefore, overall a periodic 
superfluid is a non-Landau superfluid. 

Note that a more detailed version of the above discus- 
sion for the superfluidity of a periodic superfluid can be 
found in Ref. |38j . However, an error was made in Ref. 
[38j : a non-zero critical velocity for the scenarios shown 
in Fig. [5ja,b,d) was predicted. 

The supcrfluid with SOC studied in this Letter is clearly 
also a non-Landau superfluid but for a very different rea- 
son, the lack of invariance under the Galilean transforma- 
tion. 
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